Cooper pairing reexamined 
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Abstract 

When both two-electron and two-hole Cooper-pairing are treated on an equal footing in the ladder approxima- 
tion to the Bethe-Salpeter (BS) equation, the zero-total-momentum Cooper-pair energy is found to have two real 
solutions £q S = ±2I%ujd / V 'e 2 / A + 1 which coincide with the zero-temperature BCS energy gap A = Hu>d/ sinh(l/A) 
in the weak coupling limit. Here, Hlud is the Debye energy and A > the BCS model interaction coupling pa- 
rameter. The interpretation of the BCS energy gap as the binding energy of a Cooper-pair is often claimed in 
the literature but, to our knowledge, never substantiated even in weak-coupling as we find here. In addition, we 
confirm the two purely-imaginary solutions assumed since at least the late 1950s as the only solutions, namely, 
So 3 = ±i2fiw D /Ve 2 A-l. 

The bound-state, two-particle Bethe-Salpeter (BS) 1] wavefunction in the ladder approximation, with both particle- 
and hole-propagation, for the ideal Fermi gas (IFG)-based generalized Cooper pair (CP) problem [3] is 

V(k,£) = \ G (K/2 + k,£ K /2 + E) G (K/2 - k, £ K /2 - E) x 
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Here L d is the "volume" of the d-dimensional system; K = ki + k2 is the total or center-of-mass momentum (CMM) 
and k = ^(^ — k2) the relative momentum wavevectors of the two-particle bound state whose wavefunction is ifj(k,E); 
u(|k — k'l) is the Fourier transform of the interparticle interaction, £k = E\ + Ei is the energy of this bound state 
while E = E\ — E2, and Go (K/2 + k, £/2 + E) is the bare one-fermion Green's function given by ((2], p. 72) 

y ' iX-Ei+e^- E F -ze -E x + e kl - E F + is / K> 

where = h 2 kf/2m and 9(x) is the step function, so that the first term refers, e.g., to electrons and the second to 
holes. The latter are also fermions but of positive charge +e. 

Consider first the case where holes are ignored, i.e., neglect the second term in ([2]). Note that the energy dependence 
in {l} derives from the Green's function only and therefore allows defining a new function ip^ by first writing 

1&(k,£0 = G (K/2 + k, £ K /2 + E) G Q (K/2 - k, £ K /2 - E) ^ 

which upon substitution in ([1]) yields 
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dE , 0(k' ~ M 8{k' - k F ) 

-£ K /2 -E' + e K /2+k' - E F - is -£ K /2 + E> + e K/2 _ k , - Ep - is ' ( ' 
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The energy- E' integration then leaves 

i ^ £K/2+k' + £K/2-k' + - 

which may be recognized as the Bethe-Goldstone (BG) equation [J; see Fig. 1 below. 
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Figure 1: The Bethe-Goldstone equation (left) considers only 2e-CPs. The more general BS (right) includes 2h-CPs. 

For the ideal-Fermi-gas-sea-based scenario for CPs when holes are not neglected, we assume the BCS model 
interaction u( k — k'l) = — V9Bcs( e k)&BCs{ e k') where V > is the strength of the net attraction between pair 
partners and the unit step functions Obcs{^) restrict particle or hole energies e^, £fc 2 to an energy interval of width 
2hujo around the Fermi level Ep, namely Ep < ek 1 ,£k 2 < Ep + hujp> (for particles) and Ep — hujo < ek l; ek 2 < Ep 
(for holes). Integration over energies in ^ can then be evaluated directly in the complex -E'-plane resulting in the 
following equation for the wavefunction ip-^ with zero CMM K = that is 

^ = 7? ^ ' {2z»-*Ep-e f« ~ 7? 5 "(2e k ,-2E F -£ f k ' (5) 

where £q is the K = eigenvalue energy. The single prime over the first (2e-CP) summation term denotes the 
restriction Ep < < Ep + fkupi while the double prime in the last (2h-CP) term means Ep — Hlod < £k' < Ep. 
The ordinary Cooper (or BG) problem is compared in Fig. 1 with the BS problem where electron-hole symmetry is 
restored through inclusion of 2h-CPs, represented by the second term of J5J, in addition to the 2e-CPs. Ignoring the 
second term of ([5]) gives the well-known solution [5] 

£ C = -2huj D /(e 2/x - 1) ^> -2ftw D exp(-2/A) (6) 

corresponding to a negative-energy, stationary-state bound 2e-CP, where A = VN(Ep) > with N{Ep) the electronic 
density of states for one spin. Note that a 2e-CP state for general CMM wavevector K of energy £^ is characterized 
only by a definite K but not definite relative-momentum wavevector. This alone implies that ordinary, as well as the 
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generalized CPs to be considered below, obey Bose statistics [6l [7] . Without the first summation term in ([5]) the same 
expression for the Sq of 2e-CPs follows for 2h-CPs, apart from an overall sign change. 

Since (J5J) includes 2h-CPs along with 2e-CPs, eliminating the ip^'s from leads to the BS eigenvalue equation 
for the pair energy So 

Ef+Hwu Ef 

f de f de 



e-E F - £ /2 J e — Ep — S /2 ' 

Ef Ef — Iilud 

This is precisely Eq. (7-7) of Ref. [5j (where all energies are measured from the Fermi level) and in slightly different 
form than Eq. (33.2) of Ref. [9] before assuming that So is pure imaginary. If we assume that £ can be real and < 
this eigenvalue refers to the 2e-CP "sector" in the BCS model interaction. In this case, the denominator in the first 
term never vanishes since the pole lies below Ep. A similar argument holds on considering the 2h-CP sector when 
So > and now ignoring the first integral. 

However, if we consider both particles and holes simultaneously as in ([?])■ then we must take into account that, if 
there exists a real binding energy, there is a pole in one or the other integration intervals depending on the sign of So- 
li one assumes that So/2 — ia, as in [8] or in [9] with a real both integrals are now free of singularities and can be 
integrated directly to give 
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This yields the well-known pair of purely-imaginary roots reported in Refs. [8l [9l [TO], namely 



(8) 



S ( f s = ±i2hw D / Vexp(2/ 'A)-l — * ±i2hbJr>exp(—l/X) = ±i lim A (9) 

where 

A = hu D /smh(l/\) (10) 

is the zero-temperature BCS energy gap [llj . 

A more general solution So can be obtained by assuming, without loss of generality, that it be complex, namely 

So = rexpi^ (11) 

in ([7]). Then, direct integration in both terms of ([7]) yields 

s 2 - (2tuj D y 

-{) 

Defining 



-=ln— \ 2 ' = ln[p exp iO] = \np+ id. (12) 

A Sr. 



(2hto D ) 2 /r 2 = /3 2 > (13) 
and equating real and imaginary parts of (fT2"|) leads to 

a) p= v/l + /3 4 - 2/3 2 cos20 and b) 6 = tan" 1 -J^^— . (14) 

1 — p z cos 20 

Clearly, (JT2J) means that 

2/A = lnp and = 9. (15) 
The last identity substituted in (fT4b ) implies that 1 sin 2<f> = which in turn is satisfied for 

4> = 0, ±tt/2, ±tt, ±3tt/2, • • • . (16) 



Consider the solutions with = ±7r/2 so that (TJjEh.) becomes p = \Jl + /3 4 + 2/3 2 = 1 + f3 2 > 0. Thus, the first 
relation in (fTS"]) exponentiated gives exp(2/A) = 1 + 1 and again recalling the definition (|13p for (3 as well as (TTT|) 
leads to exp(— i(/))So = r = 2fruip,l -\Z cx p(2/A) — 1 whereupon inserting <fi = ±n/2 gives precisely ©. 



3 



However, if we take the first solution (f> = of (| 16[) in (fT4a ) this becomes p = a/(1 — /3 2 ) 2 = 1 1 — /3 2 1 . Since now 
£ o = r and recalling (TT3|) one obtains the reai energy eigenvalues 



BS . 



± 2hhJD = — > ±2fhj D exp(-l/A) = ± lim A. (17) 



Vexp(2/A) + 1 



In magnitude, this is always much larger than ([6]) since e -1 / A 3> e~ 2 / x as A — » 0. In Fig. 2 we plot the ratio £q S /£q 
of the exact (07]) to the exact ([6]) (upper curve); the ratio A/£q S of exact BCS gap A (fT0|) to real BS binding energy 
£(f s found here (lower full curve); dashed curve is ratio of exact A to its weak-coupling limit, rhs of (flQ|) . Even for 
A as large as 1/2 (the Migdal upper limit [T^], marked as thin vertical line; see also Ref. [T3] p. 204.) this ratio is 
only about 1.019. Successive values ±7T, ±37r/2, • • • for (j> in (fTo) give nothing new. Solutions © are the well-known 
imaginary roots of the BS integral equation in the ladder approximation; they have been reported in Refs. [8l [9l [TO]. 
The purely real energies ()17[) found here appear to be new. 
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Figure 2: Bottom: Ratio of exact BCS gap A (|10[) to real BS binding energy found here (lower full curve). Dashed 
curve is ratio of exact A to its weak-coupling limit, extreme rhs of (|10p . Upper full curve is ratio of real exact £q S 
(fl~7|) to exact £§ ([6]). Horizontal thin dashed line marks drastic scale change. Vertical thin line is the Migdal upper 
limit pj] on A. 



In summary, we find real solutions for the binding energy of a "generalized" Cooper pair when the underlying 
BCS-type interaction is allowed to act between pairs of holes as well as of electrons through a Bethe-Salpeter equation 
that restores particle-hole symmetry around the Fermi level. The magnitude of the binding energy coincides with the 
BCS energy gap in the weak coupling regime. Finally, we note that the correct physical CP binding energies ±2A 
instead of ±A, follow when the ideal- Fermi-gas sea is replaced by a BCS-correlated sea, as in Ref. 0] Eqs. (12) and 
(13). 
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